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Abstract. A generic property of a first-order phase transition in equilibrium, and in the 
limit of large entropy per unit of conserved charge, is the smallness of the isentropic 
speed of sound in the "mixed phase". A specific prediction is that this should lead to 
a non-isotropic momentum distribution of nucleons in the reaction plane (for energies 
~ 40A GeV in our model calculation). On the other hand, we show that from present 
effective theories for low-energy QCD one does not expect the thermal transition rate 
between various states of the effective potential to be much larger than the expansion 
rate, questioning the applicability of the idealized Maxwell/Gibbs construction. Ex- 
perimental data could soon provide essential information on the dynamics of the phase 
transition. 



1. Introduction 

Heavy-ion collisions at high energy are devoted mainly to the study of strongly interact- 
ing matter at high temperature and density. In particular, a major focus is the search for 
observables related to restoration of chiral symmetry and deconfinement, as predicted by 
lattice QCD to occur at high temperature, and in thermodynamical equilibrium [ [j]]. Also, 
various effective models for low-energy QCD predict chiral symmetry restoration at high 
net baryon density and low temperature [ §]. 

Several observables for the transition to deconfined matter have been proposed, such as 
electromagnetic radiation, strangeness enhancement and equilibration, charmonium disso- 
ciation, and "irregularities" in the hydrodynamic flow pattern [ |J]. The latter, in particular, 
represents an observable that is related to the effective potential at finite temperature and 
density, because hydrodynamical expansion is driven by pressure gradients, and the pres- 
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sure is minus the value of the effective potential at the global minimum. Therefore, if 
indeed a new minimum of the effective potential opens up at some temperature or density, 
the pressure gradients should change and allow observation of that new state. Here, we dis- 
cuss the differential distribution of the momenta of nucleons in the reaction plane (spanned 
by the beam axis and the impact parameter axis) as an observable for the emergence of a 
second minimum of the effective potential, which defines a new thermodynamical state of 
hot matter. Finally, we shall discuss if a phase transition close to equilibrium is likely to be a 
reasonable approximation, how one can improve on that model, and possible consequences 
for observables. 

2. Effective Potential for the Chiral Phase Transition 

For illustration, consider the following Lagrangian for the approximately 0(4) symmetric 
chiral field 4> = (a,7t) coupled to (constituent) quarks q — (u,d): 

l =q[ifd^-g{a + iy 5 T:-n)}q + ^(d fl ad^a + d p Kd fJ n)-U(a,n) . (1) 

The potential, exhibiting both spontaneously and explicitly broken chiral symmetry, is 

X 2 

U{o,%) = —(o 2 +n 2 -v 2 ) 2 -H q c . (2) 

The vacuum expectation values of the condensates are (a) =f n and (k) — 0, where f % = 
93 MeV is the pion decay constant. The explicit symmetry breaking term is due to the finite 
current-quark masses and is determined by the PCAC relation which gives H q = f^rtS^, 
where m x — 138 MeV is the pion mass. This leads to v 2 = f 2 —m^/X 2 . The value of A, 2 = 20 
(used throughout this paper) leads to a a-mass, m 2 = 2X 2 f 2 + m 2 , equal to 600 MeV. 
For g > 0, the finite-temperature one-loop effective potential also includes the following 
contribution from the quarks: 

V,(*) = dq T J ^\og(\+e- E ' T ) . (3) 

Here, d q = 24 denotes the color-spin-isospin-baryon charge degeneracy of the quarks. 
V q (<!>) depends on the order parameter field <t> through the effective mass of the quarks en- 
tering the expression for the energy E = \/k 2 + g 2 <t> 2 . The quarks constitute the heat bath 
in which the long-wavelength modes of the chiral field, i.e. the order parameter, evolve, re- 
sulting in the finite-temperature effective potential V e g = U + V q . For rather large values of 
g, the theory exhibits a first-order (chiral symmetry restoring) phase transition. For exam- 
ple, for g = 5.5 the critical temperature is T c = 123.7 MeV at vanishing baryon-chemical 
potential /j. At the phase transition line in the T — /j plane, V e g exhibits two degenerate 
minima labeled (a) and (b) in Fig. [IJ respectively. The first is the thermodynamical state 
corresponding to restored chiral symmetry, while the second corresponds to broken sym- 
metry and smoothly approaches the physical vacuum at T = 0. The two states are separated 
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Fig. 1. An example for a finite-temperature effective potential exhibiting two states (a) and 
(b), respectively, as a function of the order parameter field <t> in a-direction and temperature 
T, for vanishing net baryon charge. See text and ref. [ Q| for details. 

by a barrier which becomes smaller as the temperature decreases, ending in a point of in- 
flection at T sp w 108 MeV (for our set of parameters and fj = 0), the so-called spinodal 
instability. At this temperature, there is no more barrier to the true ground state and the 
order parameter "rolls" down towards the state (b), a process called spinodal decomposi- 
tion. In thermodynamical equilibrium, the order parameter is localized in either one of the 
two minima, whichever is at lower energy. At the phase boundary, where the two states are 
degenerate, the system is in the "mixed phase", the expectation value of <t> being peaked 
both about <£>( a ) as well as A Maxwell construction can be performed two obtain the 
relative probabilities for the system to be in either of the two states. 

If one is interested in the hydrodynamical evolution of a system in local equilibrium, 
and described by the effective potential V e g, all one needs to know is the value of V e g at its 
global minimum, which is the value of the grand canonical potential (density) at the given 
temperature and chemical potential, i.e. minus the pressure. Once the function p{T,fu) is 
known, the evolution (in local equilibrium) of some initial condition is uniquely determined 
by the equations of ideal relativistic hydrodynamics. 



3. Hydrodynamics of Strongly Interacting Matter 

Hydrodynamics is defined by (local) energy-momentum and net charge conservation [ 



a^ v = o , 3^ = 



(4) 
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denotes the energy-momentum tensor, and M 1 the net four-current of the zth con- 
served charge. We will explicitly consider only one such conserved charge, namely the 
net baryon number. We implicitly assume that all other charges which are conserved on 
strong-interaction time scales vanish locally. The corresponding four-currents are therefore 
identically zero, cf. eq. (Q), and the conservation equations are trivial. 

For ideal fluids, the energy-momentum tensor and the net baryon current assume the 
simple form [ ||] 

r" v = {z+p)u^u v - P g" v , N^ = puf , (5) 

where e, p, p are energy density, pressure, and net baryon density in the local rest frame of 
the fluid, which is defined by Ng = (p,0). = diag(+, — , — , — ) is the metric tensor, and 
(/' = y( 1 , v) the four- velocity of the fluid (v is the three-velocity and y = (1 — v 2 ) ~ l l 2 the 
Lorentz factor). The system of partial differential equations (Q) is closed by choosing an 
equation of state (EoS) in the form p = p(e, p). 

The EoS employed in this section exhibits a first order phase transition to a Quark- 
Gluon Plasma (QGP). The hadronic phase consists of nucleons interacting via relativistic 
scalar and vector fields [ §]; the pressure is obtained from the one-loop effective potential 
for the nucleons, similar to eq. (|3]), plus that of free thermal pions. The QGP phase is de- 
scribed within the framework of the MIT-Bag model as an ideal gas of u and d quarks 
and gluons, with a bag parameter B l l A = 235 MeV, resulting in a critical temperature 
T c ~ 170 MeV at p = 0, while the critical baryon-chemical potential is p c = 1.8 GeV 
at T = 0. The first order phase transition is constructed via Gibbs' conditions of phase 
equilibrium. Thus, by construction the two states of the effective potential are coexisting 
in equilibrium: as the fluid expands, an increasing fraction of matter is transfered from 
the high-temperature state (a) to the low-temperature state (b). Supposedly, that happens 
through nucleation of bubbles of hadronic matter within the QGP matter [ ^ If that 
process occurs arbitrarily close to equilibrium, the two thermodynamical states remain de- 
generate throughout the phase transition region, such that the pressure is the same in both 
states. 



3.1. The effect of a first- order equilibrium phase transition on directed flow 

The equilibrium (first-order) transition discussed above has interesting implications regard- 
ing the hydrodynamical expansion pattern of the hot and dense matter [ ^| . Expansion of 
a perfect fluid conserves the entropy current s M = su M , as follows from eqs. (Q) upon con- 
traction with the contravariant four-velocity u v , and from the thermodynamical identities 
de = Tds + pdp, e = Ts — p + /jp. Thus, the entropy density current is proportional to 
the baryon density current iV#; each fluid element traces a path of s/p = const, in the phase 
diagram. The pressure gradient among neighbouring fluid elements, which is responsible 
for the acceleration in a given spatial direction, is given by Vp = c 2 Ve, where 



2 dp 



(6) 

.v/p=const. 



denotes the isentropic speed of sound. Fig. ^Jshows the pressure as a function of the energy 
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Fig. 2. The pressure as a function of the energy density at various 'values for the specific 
entropy (for the EoS described in the text). 

density at various values for the specific entropy s/p. (The "wiggles" in the curves are due 
to the plotting procedure and do not reflect the numerical accuracy, which is much better.) 
One clearly observes a rapid increase of p(e) within the hadronic phase, followed by a less 
rapid increase in the coexistence ("mixed") phase, and finally the transition to pure quark- 
gluon matter at very high energy density. In particular, the hotter the fluid, that is the larger 
s/p, the more does p(e) flatten out in the mixed phase. For s/p -nx>we obtain p = const, 
inbetween the two pure phases. That behavior can be understood immediately as following 
from Gibbs' conditions of phase equilibrium, s/p — » °° implies p/T — » 0, and thus in this 
limit the temperature is the only remaining intensive thermodynamic al variable the effective 
potential (and the pressure) can depend on. Now, in phase equilibrium T = T C = const, and 
so the pressure p(T c ) = p c is constant as well. In other words, the isentropic speed of 
sound c 2 s must be small during the time when a fluid element's trajectory through the phase 
diagram coincides with the phase coexistence line (in the T — p plane). 

If that prediction based on the equilibrium phase diagram is indeed relevant for colli- 
sions of heavy ions, it should be visible in the excitation function of the so-called directed 
in-plane flow [ H ^|, (p x /N). That observable is defined as the average momentum per 
nucleon in impact parameter (jt— ) direction [ |oJ, and due to the kinematics of the collision 
is particularly sensitive to the early stage of the reaction, where the phase transition might 
occur. A more refined observable is the distribution of momentum in the reaction plane, 
which is proportional to the triple-differential cross section for the process AB — > p + X, 
da/d 3 /?. To zeroth approximation, it can be obtained by integrating the baryon current 
over the freeze-out hypersurface of the nucleons, E'': 




(7) 



In what follows we assumed for simplicity that freeze-out, or at least onset of strong dis- 
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sipative correction to the perfect-fluid energy-momentum tensor, occurs on an equal-time 
hypersurface in the center of mass frame, dE A ' = (d 3 x,0). Furthermore, we integrated the 
differential cross section over the momentum component p y , perpendicular to both the im- 
pact parameter vector and the beam axis. 
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Fig. 3. Distribution of the time-like component of the net-baryon four-current in momen- 
tum space (p x — p z plane). Pb+Pb collisions at b = 3 fm, E^'" b = 8A GeV (left frame) and 
40A GeV (right frame), respectively. 

Typical distributions of nucleons in momentum space are depicted in Fig. ||] [ pLl| ]. The 
left frame corresponds to relatively low collision energy and specific entropy s/p w 10. 
Therefore, cl is not small in that case even if the energy density at the center exceeds 
~ 0.5 GeV/fm 3 , and the state (a) of the effective potential should be populated according 
to the equilibrium EoS (see Fig fy. So, since c 2 is not small, (energy-) density gradients in 
coordinate space reflect in non-vanishing pressure gradients Vp, which in turn act to make 
the distribution of nucleons in momentum space more or less isotropic. The net momentum 
in x-direction at any fixed longitudinal momentum p z is obviously rather small, as positive 
and negative contributions largely cancel [ [l2|] . 

On the other hand, Fig. ^ also shows that the structure of the flow changes qualitatively 
at higher energy, E^ b w 40A GeV, where on average s/p « 20. The fact that c 2 s is now 
smaller than at lower energy, by a factor of two in the present model, actually prevents more 
isotropic redistribution of the baryon number in momentum-space [ |TJ. The distribution 
is clearly very different from that at the lower energy, with almost no nucleons in the upper 
left or bottom right quadrants where p x ■ p z < [ [y]] . 

4. The nucleation rate, and supercooling down to the spinodal 

The Gibbs/Maxwell construction relies on the assumption that the thermal barrier pene- 
tration rate is much larger than the local expansion rate. In a slowly expanding system, 
the phase transition would proceed through the nucleation of bubbles of the "true vacuum" 
state via thermal activation [ Hj . The nucleation rate per unit volume per unit time is 
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expressed as 

r = fP e - Fl 'l T , (8) 

where Ff, is the free energy of a critical bubble and where the prefactor T provides a mea- 
sure of the saddle point of the Euclidean action in functional space. Our analysis will 
concentrate on the exponential barrier penetration factor, and we will approximate ¥ by 
T 4 , which is obtained from dimensional considerations at ji — 0. In order to determine the 
role of bubble nucleation in the evolving system and to be able to compute the decay rate 
r, it is necessary to study the critical bubble and some of its features. The critical bubble 
is a radially symmetric, static solution of the Euler-Lagrange field equations that satisfies 
the boundary condition <t>(r — > °°) — > ^(a)- Energetically, this boundary condition cor- 
responds to an exact balance between volume and surface contributions which defines the 
critical radius R = R c . The critical bubble is unstable with respect to small changes of its 
radius. For R < R c , the surface energy dominates, and the bubble shrinks into the false vac- 
uum. For R > R c , the volume energy dominates, and the bubble grows driving the decay 
process. 

The critical bubble can be found by minimizing the free energy 



F b (®,T) =4JC / r z dr 



1 /d<i> N 2 



2 ld7» +V ^ T) 



(9) 



with respect to the field <t>. This can be performed numerically without further approxima- 
tion [ [|]. Fig. |] shows the resulting decay rate, divided by an assumed expansion rate of 
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Fig. 4. r 1//4 divided by an assumed expansion rate of H — 1 fm~' as a function of temper- 
ature (left frame). Correlation length of the chiral order parameter in state (a) and in state 
(b). The latter corresponds to an equilibrium transition. 

H = 1 fm _1 , which represents a rough estimate for the expansion rate of the three-volume 
of a comoving fluid element, dV = d£ • u, in high-energy collisions. Y X I A specifies the 
inverse spatial and temporal scale of thermal field fluctuations into the broken symmetry 
state. Clearly, for the above effective field theory T 1 / 4 is not much larger than H, even if 
one takes into account that both factors are only known approximately. While the results of 
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Fig. |4| do not provide an unambiguous prediction for the fate of the supercooled state, they 
do make it appear likely that the rapidly expanding system has an appreciable probability 
of remaining in the restored symmetry phase even close to the spinodal instability. Thus, at 
least some fraction of all heavy ion events should show traces of this non-equilibrium tran- 
sition. The familiar idealized Gibbs/Maxwell construction of equilibrium thermodynamics 
may not be appropriate for the description of phase transitions in high-energy heavy ion col- 
lisions. If so, the anisotropy of the differential nucleon cross section in the reaction plane, 
predicted in the previous section on the basis of the equilibrium phase diagram, should be 
absent. Rather, the small net momentum in impact parameter direction (around p z ~ 0), as 
observed experimentally [|l3|] at the top BNL-AGS energy E^" b w 10A GeV, should remain 
essentially zero or diminish even further. That is because the phase equilibrium with small 
velocity of sound c 2 s does not occur in a rapid out of equilibrium transition, where the order 
parameter remains localized in state (a) of the effective potential down to the spinodal. 

Instead, other experimental observables of the phase transition may emerge. Consider 
the behavior of the correlation length of the <t>-field. It is obvious from Fig. [j] that the 
curvatures of the effective potential at the minima <!>(„) and <&(u are rather different for 
T < T c . Thus, the effective mass, m^ ff = d 2 V e ff/d<I> 2 , and the correlation length, I = l/m e g, 
of the field will also be different. In particular, if the rapidly expanding system supercools 
appreciably and if the order parameter is trapped in the metastable state, one can expect a 
clear increase of /. The right frame of Fig. Q shows I in the states (a) and (b), respectively, 
as a function of the degree of supercooling. The correlation length at (a) increases as that 
minimum disappears. By contrast, there is a smooth decrease in the correlation length at 
the true global minimum, corresponding to an equilibrium transition. In the present model 
h a \ can exceed Ir^ by as much as a factor of 2 — 3, depending on the degree of supercooling. 
The precise values of / depend on the specific effective model adopted, but the qualitative 
observation that l^s increases as the system approaches the spinodal instability is general. 

5. Conclusions 

We discussed generic properties of the effective potential for a first-order phase transition 
in equilibrium, in particular the vanishing of the isentropic velocity of sound in the limit of 
very large entropy per unit of conserved charge (baryon charge in our case). As an example 
for how such properties of the effective potential show up in the dynamics of hot and 
dense QCD matter we discussed the so-called directed flow of nucleons, i.e. the differential 
cross section in the reaction plane. Furthermore, we showed that the forthcoming results 
of the Pb+Pb reactions at E^ b = 40A GeV can test the applicability of the picture of hot 
QCD-matter as a heat-bath with small isentropic speed of sound (mixed phase) to heavy- 
ion collisions in that energy domain. If it holds true, our model calculations predict an 
increase of the directed net in-plane momentum as compared to top BNL-AGS energy, 
and a nonisotropic momentum distribution around midrapidity. During the assumed first- 
order equilibrium phase transition pressure gradients along isentropes are too small to work 
towards a more isotropic momentum distribution. 

The experimental test of the equilibrium phase transition picture would also provide 
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valuable constraints on effective field theories for low-energy QCD. As mentioned above, 
for the present model thermal barrier penetration from state (a) to state (b), and vice versa, 
is not materially larger than rough estimates for the expected expansion rates, thus leading 
to the possibility of significant supercooling, close to the spinodal instability. 

The experimental observation of supercooling effects and spinodal decomposition would 
also be important as a matter of principle. Ideally, signatures of phase transitions should 
be order parameter related and should reveal properties of the equilibrium phase diagram. 
While the spinodal instability is not part of the equilibrium phase diagram, it is rather close. 
Familiar experiments in condensed matter physics on a wealth of hysteresis phenomena (i.e. 
the analogue of superheating and supercooling) make it clear that spinodal instabilities can 
be studied on "macroscopic" time scales [ [14]]. The spinodal instability, if found, would 
certainly be among the most direct information which relativistic heavy ion collisions can 
provide regarding the QCD phase transition. 
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